Spin Degree of Freedom in a Two-Dimensional Electron Liquid 
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We have investigated correlation between spin polarization and magnetotransport in a high mobility 
silicon inversion layer which shows the metal-insulator transition. Increase in the resistivity in a 
parallel magnetic field reaches saturation at the critical field for the full polarization evaluated from 
an analysis of low-field Shubnikov-de Haas oscillations. By rotating the sample at various total 
strength of the magnetic field, we found that the normal component of the magnetic field at minima 
in the diagonal resistivity increases linearly with the concentration of "spin-up" electrons. 
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A metal-insulator transition (MIT) observed in Si 
metal-oxide-semiconductor field-effect transistors 
(Si-MOSFET's) and other systems §-§] attracts a great 
deal of attention since it seems to contradict an impor- 
tant result of the scaling theory by Abrahams et al. 0] 
that the conductance of a disordered two-dimensional 
(2D) system at zero magnetic field goes to zero for T — > 0. 
In the metallic phase in Si-MOSFET's with high peak 
electron mobilities of ^ pca k ^ 2 m 2 /V s, the diagonal 
resistivity p xx shows a sharp drop with decreasing tem- 
perature from about 2 K JO. Recent experiments 
show that magnetic fields applied parallel to the 2D plane 
suppress the low temperature metallic conduction in Si- 
MOSFET's. Since the parallel magnetic field does not 
couple the orbital motion of electrons, this fact suggests 
an important role of the spin of electrons. However, the 
mechanism of the conduction in the anomalous metallic 
phase is not clear yet. 

The 2D systems that show the MIT |]-|| are char- 
acterized by strong Coulomb interaction between elec- 
trons. The mean Coulomb energy per electron U = 
(7rA r s ) 1 / 2 e 2 /47T£o/t is larger than the mean kinetic energy 
K = nh 2 N s /m* by an order of the magnitude around 
the critical point for the MIT. Here, N s is the electron 
concentration, n is the relative dielectric constant at the 
interface, and m* is the effective mass of electron. It 
is estimated that U — 120 K, K = 14 K and the ra- 
tio r s = U/K = 8.3 for k = 7.7 and m* = 0.19m„ at 
N. 



1 x 10 15 m- 2 in Si-MOSFET's. The ground state 
of the insulating phase of high mobility Si-MOSFET's 
is considered to be a pinned Wigner solid (WS) @,|n). 
Magnetic field dependence of the thermal activation en- 
ergy observed for various angles of the magnetic field was 
essentially explained by a model based on magnetic inter- 
actions in the pinned WS fl2|;[l3| . Although the quantum 
fluctuations change the 2D system into a liquid at higher- 
N s , electron-electron (e-e) interaction is expected to be 
still important. 

In the conduction band of silicon, the spin-orbit in- 
teraction is negligible and the spin polarization p = 
(JVf — Ni)/N s can always be given in the direction to 



the magnetic field. Here JVj and N± are the concentra- 
tions of electrons having an up spin and a down spin, 
respectively (N s = iVf + iV|J. In the present work, we 
investigate the low temperature conduction in a high mo- 
bility Si-MOSFET for various values of p. 

The Si-MOSFET sample used has a peak electron mo- 
bility of /i pcak = 2.4 m 2 /V s at N s = 4 x 10 15 m~ 2 and 
T = 0.3 K. It has a Hall-bar geometry of total length 
3 mm and width 0.3 mm. The estimated Si02 layer 
thickness is 98 nm. The diagonal resistivity p xx as well 
as the Hall resistivity p xy was obtained from the linear 
relationship in the I- V characteristics using a four-probe 
dc method. The source-drain current and the potential 
difference across the two potential probes separated by 
1.5 mm were limited below 10 nA and 0.4 mV, respec- 
tively. The sample was mounted on a rotatory thermal 
stage in a dilution refrigerator together with a GaAs Hall 
generator and a Ru02 resistance thermometer calibrated 
in magnetic fields |l3[ . 

As shown in Fig. Ill, the MIT is clearly observed in a 
zero magnetic field. The critical value of p xx at the MIT 
is 55 k57. It is close to p c ~ 2h/e 2 reported by Kravchenko 
et al. [EJ. A high magnetic field of 9 T parallel to the 2D 
plane drastically increases p xx and eliminates the positive 
temperature dependence even in the low-p xx region jli} . 
The disappearance of the MIT due to the parallel mag- 
netic field seems to resemble that induced by increasing 
disorder reported by Popovic et al. 0] . 

In non-interacting degenerate 2D Fermi gases, the spin 
polarization p increases linearly with the total strength 
-Btot of the magnetic field for p < 1 and reaches p = 1 at 
the critical magnetic field of B c = 2ttH 2 N s / \iB9v9*fn* '■ 
Here, pb is the Bohr magneton (= he/2m e ). The valley 
degeneracy g v is 2 on the (001) surface of silicon and the 
effective g-factor g* is 2.0 in the conduction band in sili- 
con jl5| . According to Landau's Fermi liquid theory, the 
e-e interaction changes the effective g- factor and mass 
of electrons (or quasiparticles) into enhanced values de- 
noted by <7fl and tofl, and reduces the critical magnetic 
field. 

The product of 5fl and m FL in (001) Si-MOSFET's 
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FIG. 1. Diagonal resistivity as a function of electron con- 
centration for different temperatures in a zero magnetic field 
(closed symbols) and in a parallel magnetic field of 9 T (open 
symbols). 
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FIG. 3. The enhancement factor a = <?FL?riFL/0.38m e vs 
the strength of the electron-electron interaction. Inset shows 
the critical magnetic field B c calculated using Eq. (hi). 



can be determined from the Shubnikov-de Haas oscilla- 
tions in tilted magnetic fields p6[ . Figure || shows the 
diagonal conductivity a xx = p xx /{p xx 2 + Pxy 2 ) divided 
by eN s for various angles of the magnetic field to the 2D 
plane. The normal component B± of the magnetic field 
is fixed at 0.80 T. The upper abscissa is given by a half of 
the Landau level filling factor v = N s cj)o/B± (</>n = h/e) 
taking into account the valley degeneracy g v = 2. We 
ignore the valley splitting energy since no indication is 
observed at v — odd. For B to t = 0.80 T, a xx /eN s takes 
minima at even- numbers of v/2 where the Fermi energy 
is located at the middle of a gap between Landau levels 
having different Landau quantum number and spin index. 
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FIG. 2. Shubnikov-de Haas oscillations at B± 
and T — 0.31 K for various Btot- 
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This energy gap AE even = heB±/m F i J - gFLl^BB to t de- 
creases with Btot and the minima in a xx /eN s at even- 
numbers of is/2 fade away at higher- B to t due to the 
level broadening. On the other hand, the spin splitting 
AE dd = SFL^B^tot increases with Btot and the minima 
in <r xx /eN s at odd-numbers of v/2 become clearer. We 
have evaluated the critical angle 8 C of the magnetic field 
to the 2D plane where AE a dd becomes equal to AE even 
from comparing the depth of the minimum in a xx /eN s at 
an odd- number oiv/2 with the average of those at v/2 — 1 
and v/2 + 1 for various values of B to t- We assumed that 
the level broadening is a smooth function or independent 
of v. Thus we have <7fl?tifl = m eB±/ B to t = m e sm6 c . 

In Fig. ||, we show the ratio of <7fl?tifl to g*m* in 
the non-interacting system. The enhancement factor 
a = gFLTOFL/0.38m e increases monotonically with the 
strength of the e-e interaction. Results on much lower 
mobility samples measured in higher magnetic fields by 
Fang and Stiles |ll| are also plotted. No B^-dependence 
is found between our results at 0.8 T and at 1.2 T. We 
believe that these values of B± are small enough and 
gpLmpL was determined in the limit of B± = 0. While 
B± may cause a small oscillation of the effective g-factor 
with the maxima gp^ x at v/2 = odd and the minima 
3fl N at = even ]17| ], 3fl"1fl determined from the 
critical angle for AE a dd = AE even is related to the av- 
erage of gpL^ and <?fl N . The present method cannot 
be used in the large-r s region where gFxjn^/me exceeds 
unity (a > 2.63). We use a simple extrapolation function 

a - 1 = 0.2212r s + 0.003973r s 2 . (1) 

represented by the solid line in Fig. ^|. The critical mag- 
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FIG. 4. Resistivity vs B tot /B c at T = 0.21 K and B ± = 
for various N s . Dashed line represents the critical value of 
Pxx at the MIT determined from temperature dependence be- 
tween 0.21 K and 0.91 K. 

netic field B c = 2-kTi 2 'N s / /LtB^gFLm-FL is calculated using 
(||) as shown in the inset to Fig. [| 

In Fig. ||, p xx at T — 0.21 K in the parallel magnetic 
field (B± = 0) is shown. The diagonal resistivity p xx 
increases with B to t in the low-B t ot region, but it takes 
almost constant values in the high-B to t region. The sat- 
uration of p xx occurs around B tot — B c at which the spin 
polarization is expected to be completed. The result in- 
dicates that the mixing of the different spin states causes 
the reduction of p xx in the low-i?tot region. The dashed 
line represents the critical value of p xx at the MIT. It 
was tentatively determined from the sign of the change 
in p xx from T = 0.21 K to 0.91 K. The positive tem- 
perature dependence of p xx below 1 K was not observed 
for B to t > B c . The positive T-dependence of p xx in the 
low-i?tot region may arise from the scattering related to 
the spin degree of freedom. 

The most important finding has been obtained by ro- 
tating the sample in the magnetic field. Results at 
N s = 1.46 x 10 15 m~ 2 are shown in Fig. ||. We focus on the 
p xx oscillation that appears in the low-i?j_ region, while 
the p xx (B±) curves for different B to t approach each other 
in the higher- B± region (not shown in Fig. |^) and con- 
verge on a single curve for B±/N s ()>o > 0.5 with two deep 
minima at B±/N s (j> = 0.5 (u = 2) and B±/N s (j>o = 1.0 
(y = 1) as reported by Kravchenko et al. Unlike 
ordinary Shubnikov-de Haas oscillations that show p xx 
minima at fixed points with v = integer, the p xx minima 
indicated by the black arrows and the white arrows in 
Fig. shift toward high-Bj^ side as i? to t increases. Such 
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FIG. 5. Diagonal resistivity vs normal component of the 
magnetic field at T = 0.31 K and N s = 1.46 x 10 15 m~ 2 for 
several values of B to t- 

peculiar behavior was also observed for lower- iV s down 
to 1.02 x 10 15 m- 2 . 

In Fig. |(| the positions of the p xx minima for various 
N s are shown. Both B± and B to t are normalized by N s (j)Q 
and B c , respectively. The closed (open) symbols corre- 
spond to the black (white) arrows in Fig. | The value of 
B±/N s (f)o increases linearly with B to t/B c for B to t/B c < 
1, but saturates for Btot/B c > 1 where the spin polariza- 
tion is expected to be completed. We link this behavior 
to the concentration iVj of electrons having an up spin. 
Assuming that a change in p by B± is negligible in the 
low-i?j_ region, we have N^/N s = (1 + B tot / B c )/2 for 
B tot /B c < 1 and N^/N s = 1 for B tot /B c > 1. All the 
data are distributed along the solid line or the dashed 
line on which the number of the flux quantum per "spin- 
up" electron takes a constant value of B±/N^tf>o = 0.265 
or 0.172. 

One may attribute the observed p xx minima to the 
Shubnikov-de Haas oscillation produced by spin-up elec- 
trons. Actually, the effective Landau level filling factor 
of spin-up electrons A|0o/ B± at the p xx minima is close 
to an even number of 4 or 6. Note that we have the 
valley degeneracy g v = 2 in this system. It is expected 
that p xx also takes minima at even-numbers of the ef- 
fective Landau level filling factor of spin-down electrons. 
However, we observed no feature even for Ni<fio/B± = 2 
represented by the dotted line in Fig. ^. The disappear- 
ance of p xx minima related to spin-down electrons im- 
plies that spin-down electrons do not contribute to the 
conduction owing to localization or large scattering rate. 
If that is the case, one should accept an unusual negative 
dependence of the mobility of spin-up electrons on the 
carrier concentration iVj since the value of p xx steeply 
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FIG. 6. Positions of the p xx minima for various 
N s are plotted. Solid line and dashed line represent 
B ± /N 1( j>o = 0.265 and B^/Nrfo = 0.172, respectively. Dot- 
ted line represents B±/Ni(j>o = 0.5. 



increases with the spin polarization. The p xx minimum 
at B±/Ni<po ~ 0.265 appears even in the insulating re- 
gion with p xx ~ 10 6 fi, where the Landau level separa- 
tion is considered to be entirely smeared out by the level 
broadening. Further consideration is required to explain 
the results consistently. 

Finally, we also propose a quite different interpreta- 
tion of the p xx minima. In the previous work fl^ , |l3| , we 
have investigated the activation energy in the insulating 
phase which depends on B± and -Btot- Oscillatory be- 
havior linked with B±/N g (f>o was explained using a model 
based on the Aharonov-Bohm effect on the exchange of 
electrons localized in the pinned Wigner solid |19|. This 
model may be partly applied to the present case if one 
assumes that spin-up electrons are localized and form an 
ordered state in the real space like a Wigner solid. Then, 
the £?tot-dependence of p xx can be understood as a result 
of the change in the concentration of spin-down electrons 
which carry the charge. While we at this stage do not 
identify the origin of the p xx minima specifically, we con- 
sider that the effect of the magnetic flux on the system 
of spin-up electrons can cause an oscillatory perturba- 
tion in the mobility of spin-down electrons through some 
interactions. 

In summary, we have studied the magnetotransport in 
a high mobility Si-MOSFET at electron concentrations 
around the MIT. The relationship between the resistiv- 
ity and the spin polarization indicates that the mixing 
of the different spin states is important for the metallic 
conduction. We also found the p xx minima in the in- 
dependence, which should be related to the concentration 
of "spin-up" electrons. 
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